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We investigate three-dimensional 0(N) spin models driven with a uniform velocity over a random 
held. Within a spin-wave approximation, it is shown that in the strong driving regime the model 
with N = 2 exhibits a quasi-long-range order in which the spatial correlation function decays in 
a power-law form. Furthermore, for the cases that N = 2 and 3, we numerically demonstrate 
a non-equilibrium phase transition between the quasi-long-range order phase and the disordered 
phase, which turns out to resemble the Kosterlitz-Thouless transition in the two-dimensional pure 
XY model in equilibrium. 
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I. INTRODUCTION 

Thermodynamic systems can exhibit remarkably com¬ 
plex phase transition dynamics when they are driven by 
an external force. Well-studied examples include ferro¬ 
magnetic systems and crystalline systems that are driven 
by steady shear GL[3, oscillating external fields @,[U, and 
athermal noise [5=, 6J. The general treatment that con¬ 
cerns how the driving forces affect the nature of phase 
transitions has not been established yet. To identify and 
classify various types of phase transitions peculiar to out 
of equilibrium systems can give valuable insights for de¬ 
veloping non-equilibrium statistical mechanics. 

The ferromagnetic and crystalline phases are charac¬ 
terized by long-range order (LRO), wherein the order 
parameter takes a finite value. As a qualitatively dif¬ 
ferent behavior, the two-dimensional XY model exhibits 
quasi-long-range order (QLRO), wherein the order pa¬ 
rameter remains zero but the spatial correlation decays 
in a power-law form at low temperatures Q. Further¬ 
more, at the transition between the QLRO phase and the 
disordered phase, which is called the Kosterlitz-Thouless 
(KT) transition, there is no singularity in thermody¬ 
namic quantities such as specific heat and susceptibility 
[!j. This peculiar type of phase transition has attracted 
considerable research attention. However, the possibil¬ 
ity that a driving force causes the KT transition has not 
been clearly discussed or established. Consequently, the 
question arises as to whether there is a system that does 
not exhibit the KT transition in equilibrium, while it 
does exhibit the transition in the presence of a driving 
force. To the best of our knowledge, such a system has 
not been reported thus far. In this paper, we consider 
three-dimensional O (N) spin models driven with a uni¬ 
form velocity over a random field as examples of such sys¬ 
tems. This “non-equilibrium KT transition” may define 
a novel type of universality class, wherein the interplay 
between disorder and driving plays a crucial role. 

The dynamics of an ordered system driven by an ex- 
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ternal force over a random substrate have been a topic 
of intensive research in statistical physics. A rich va¬ 
riety of complex phenomena results from the interplay 
between elasticity, quenched disorder and driving. The 
best-known example of such systems may be driven vor¬ 
tex lattices in dirty superconductors |9l111] . In such sys¬ 
tems, impurities and crystalline defects act as random 
pinning potentials for vortex lines. There are other well- 
studied systems, e.g., charge-density waves [l2] or col¬ 
loids ra driven by an external field. The driven vor¬ 
tex lattice systems exhibit dynamical reordering tran¬ 
sition for a finite driving velocity. In equilibrium, the 
vortex lattice system has two different phases 0 - In 
the strong disorder regime, the vortex glass phase is real¬ 
ized, in which the spatial correlation function of the dis¬ 
placements of the vortices exponentially decays. In the 
weak disorder regime, the Bragg glass phase is realized, 
in which the spatial correlation function exhibits power- 
law decay [l5j| . When the vortex glass is driven with a 
finite velocity, the effective disorder that the vortices ex¬ 
perience becomes weaker because the random potential 
varies rapidly in a moving frame. At a large driving veloc¬ 
ity, a first-order phase transition to the Bragg glass phase 
from the vortex glass phase occurs [ 1(1 Hlj . While the 
vortex lattice systems that exhibit the first-order phase 
transition in the absence of the disorder are well studied, 
very few studies have focused on a driven system over 
a random substrate that exhibits a second-order phase 
transition in the absence of disorder. 

In this paper, we consider the nature of the dynami¬ 
cal reordering transition of three-dimensional O (N) spin 
models when they are driven with a uniform velocity over 
a random field. We show that the models with N = 2 
and 3 exhibit QLRO at low temperatures and that the 
transition from the QLRO phase to the disordered phase 
resembles the KT transition. This paper is organized as 
follows: In Sec. El we introduce the models and review 
their behavior in equilibrium. In Sec. ED for the case 
that N = 2 (XY model), we show that this model ex¬ 
hibits QLRO at low temperatures by using the spin-wave 
approximation. In Sec. m the phase diagram for N = 2 
and 3 with respect to temperature, disorder, and driv¬ 
ing velocity is determined by means of numerical exper- 
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iments. We calculate the specific heat and show that it 
exhibits no singularity at the transition point. In Sec.fVl 
we summarize our results. We also discuss nematic liquid 
crystals flowing in a random medium as an experimental 
realization of our model. 


II. MODEL 


Let {4> a }a =1 be an TV-component (TV > 2) real vector 
field. The Hamiltonian of the three-dimensional O (TV) 
models with a quenched random field is given by 


H[<fi: h] = J d 3 r 


\K\V(f>(r )\ 2 - 7 -\(t>(r )| 2 


+ ^|0(^)| 4 - Hr ) ■ 


(1) 


where K. r, and g represent positive parameters. The 
quenched random field h(r) obeys a Gaussian distribu¬ 
tion with ( h(r)) = 0 and {h a (r)h^(r')) = hl5 a pS(r—r'). 
The time evolution of the field <f>(r,t) is described by 


d(j> a {r , t) 
dt 


+ v- V(j> a {r,t) = -T 


SH[cf>; h] 
S(j) a (r, t) 


+ V a (r,t), (2) 


where v = ve x denotes a parameter independent of 
(r,t), and 7j a (r,t) represents thermal noise that satis¬ 
fies (r] a (r, t')) = 2 YT5 a pS{r — r')5(t — t'). This 

equation describes the relaxation dynamics of the ordered 
system that is driven with a uniform velocity over the 
quenched random field. Examples of such systems in¬ 
clude nematic liquid crystals flowing in a random medium 
such as an aerogel or a porous substrate, where the 
random field corresponds to the random anchoring that 
occurs due to the complicated surface structure of the 
porous substrate. We investigate the non-equilibrium 
steady states of the models with TV = 2 (XY model) 
and TV = 3 (Heisenberg model). We call these models as 
driven random field O(TV) models (DRFO(TV)Ms). 

Let us review the behavior of these models in equi¬ 
librium (v = 0). The random field O(TV) models 
(RFO(TV)Ms) are one of the simplest disordered systems 
in which spins with O(TV) symmetry are coupled to a 
quenched random field. The Imry-Ma argument and “di¬ 
mensional reduction” property state that LRO with con¬ 
tinuous symmetry breaking is destroyed by an infinitesi¬ 
mally weak random field below four dimensions (l3 - H^ | . 
The absence of LRO below four dimensions was also 
rigorously proved by Aizenman and Wehr [l9| . While 
the theoretical description that concerns the existence 
of LRO has been established, whether QLRO exists or 
not in three dimensions is a more subtle problem. From 
the analogy of the Bragg glass phase in a vortex lattice 
system [15(, theoretical studies based on the renormal¬ 
ization group analysis have suggested the existence of 
QLRO for the three-dimensional random field XY model 
and the random anisotropy Heisenberg model [2p| . [2lj . In 


order to confirm this predicted QLRO, numerical simula¬ 
tions have also been conducted [22|,[23|. However, defini¬ 
tive evidence for the existence of QLRO has not thus far 
been obtained because the correlation length rapidly in¬ 
creases when the strength of the random field decreases. 
Moreover, in Ref. [HI, the experimental investigation of 
nematic liquid crystals in aerogels did not lead to the 
observation of any QLRO. Recently, more sophisticated 
renormalization group studies have neg ated the existence 
of QLRO in three dimensions |25l [26| . In summary, to 
the best of our knowledge, we can conclude that three- 
dimensional RFO(TV)Ms with TV > 2 do not exhibit any 
phase transitions. 


III. SPIN-WAVE APPROXIMATION 


Let us calculate the spin correlation function for the 
XY model (TV = 2) by using the spin-wave approx¬ 
imation at zero temperature. The magnitude of the 
spin is assumed to be fixed to unity. The order pa¬ 
rameter field is represented by 4>(r) = (0 1 (r), </> 2 (r)) = 
(cos0(r), sin0(r)). The Hamiltonian is rewritten as 


i/V(W(r)) 2 - h(r )cos(0(r) - £(r)) , 

(3) 

where the random field is written as h(r) = 
(h 1 (r), h 2 (r)) = (h(r)cos£(r),h(r)sm£(r)). The equa¬ 
tion of motion at zero temperature is given by 



d0(r,t) 
dt 


+ v d9( Q^ =T[KV 2 6{r,t) 
-h(r)sin(0(r,t) - £{r))]. 


( 4 ) 


We define the Green function G(r) by its Fourier trans¬ 
form G(q) = ( TKq 2 + ivq x )~ 1 . Subsequently, the formal 
solution for the steady state is given as 


0(r) = J d r'G(r — r') 

xT {—Ti 1 (r , )sin0(r / ) + h 2 (r')cos9(r')} . (5) 

The mean square relative displacement 9(r-2)) 2 ) 

is calculated as 

<(#(»■!)-0(r a )) 2 ) 

= J dVdV' {G{ ri - r') - G(r 2 - r')} 

x {G( ri - r") - G(r 2 - r")} 
xT 2 { (h 1 (r')h 1 (r'')sm9(r')sm9(r")) 

—2(Ti 1 (r')Tj 2 (r // )sin0(r / )cos0(r / ')) 

+(h 2 (r')h 2 (r")cos9(r')cos9(r")} }. (6) 

We use factorization approximations such as 
{h a (r)hP(r')sm9(r)sin9(r')) ~ (/i“(r)/i /3 (r / )) 
x (sin0(r)sin0(r / )) [13 ■ This factorization is justified 
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when the correlation length of 6(r), which is denoted 
by £, is much larger than that of the random field £r. 
We will consider a self-consistent condition leading to 
£ £r after the calculation of the mean square relative 
displacement. If this condition is satisfied, we have 

m ri )-0(r 2 )) 2 ) 

= T 2 h 2 0 J dr' {G(n - r') - G(r 2 - r')} 2 . (7) 


Substituting the explicit form of the Green function, we 
have 


m ri )-e(r 2 )) 2 ) 

_ 9r 2,2 f d 3 g 1 — cos{q ■ (ri — r 2 )} 

°J (2tt) 3 r2/TV + v 2 g2 ‘ 1 J 

From Eq. (|8|) . we obtain the asymptotic behavior over 
a large distance r^FK/v, 


mr)-e( 0)) 2 > 


• IX 

47 tKv 

r hi 

. 2nKv 


In r, 
In r, 


(r || r), 
(r _L n). 


( 9 ) 


The detailed calculation is presented in the Appendix A. 
This logarithmic dependence on the distance is similar 
to that of the KT transition in the two-dimensional pure 
XY model [H. 

Let us consider a self-consistent condition leading 
to £ £_r. We define the correlation length £ by 

((0(£) — 6( 0)) 2 ) ~ 1. The correlation length of the 
random field £# is equal to a cut-off for short length 
scale A -1 because the correlation function is given by 
(/i“(r)h /3 (r')} = hl5 a p5(r — r'). If we suppose £ £# 

as an ansatz, we obtain from the above calculation 


£ 

A- 1 


exp 


2ttKv 


( 10 ) 


which becomes infinitely large in the weak disorder limit 
or in the strong driving limit. This argument suggests 
that £ £_r when 


2nKv 

~Thl 


> i. 


(ii) 


Thus, if the condition Eq. m is satisfied, the factoriza¬ 
tion approximation used in the calculation of Eq. 0 is 
justified. In the next section, we will confirm that this 
condition is satisfied in the region of QLRO. 

The correlation function C(r) = (4>(r) ■ </>(0)) is calcu¬ 
lated as follows: 


G(r) = ( e ‘(»(»O-0( o ») 


= exp 


o° 

E^ n ( (0 w- 0(o))ri )c 


( 12 ) 


where ((...) n ) c denotes a n-th cumulant. If one admits 
the factorization approximation noted above, it is shown 


that all higher cumulants approximately vanish because 
the random field h a obeys a Gaussian distribution. Thus, 
we have G(r) = e -1/,2 ^ e ( r )-0(o)) 2 }_ FromEq. @, we have 


G(r) oc 


r “H, (r || v), 

r~a± s ( r J_ v ^ 


(13) 


where the exponents are or|| = Thg/(8nKv) and o_l = 
r/ig/ (AttKv). Therefore, we have shown that DRFO(2)M 
exhibits anisotropic QLRO at low temperatures, in which 
the spin-wave approximation is valid. 


IV. NUMERICAL RESULTS 

We investigate the transition between the QLRO phase 
and the disordered phase by numerically solving Eq. 0. 
The calculation is implemented in a moving frame with 
velocity v. The periodic boundary conditions and free 
boundary conditions are imposed for the directions per¬ 
pendicular and parallel to the driving velocity v 1 respec¬ 
tively. The random field is continuously generated in the 
front boundary of the simulation box, and it moves with 
velocity —v. The detailed method of the numerical calcu¬ 
lation is presented in the Appendix B. Time integration is 
performed by employing the Euler method. The param¬ 
eter values are fixed as K = 1, T = 1, r = 5, and g = 10. 
We set the time and space discretization as 5t = 0.005 
and Sx = 1, respectively. 


A. XY model 

We first calculate the spin correlation function G(r) 
for the XY model (N = 2). We start from a random 
initial condition and obtain a steady state after a suffi¬ 
ciently long time. The correlation function is calculated 
from a field configuration <p(r) of the steady state. For 
a sufficiently large system, the average with respect to 
the random field and the thermal noise can be replaced 
by the spatial average because of the self-averaging prop¬ 
erty. We also take the time average, which is equivalent 
to the ensemble average with respect to the random field 
and the thermal noise. The detailed explanation for the 
method of the numerical calculation of the correlation 
function is presented in the Appendix B. The correla¬ 
tion functions for different values of temperature are dis¬ 
played in Fig. [T] The upper (a) and lower (b) panels 
depict C(r) for the directions parallel and perpendicular 
to the driving velocity v , respectively. In order to verify 
the finite size effect, we calculate C(r) for system sizes of 
60 3 , 100 3 , and 150 3 . The system size dependence of G(r) 
for the perpendicular direction is larger than that for the 
parallel direction because of the periodic boundary con¬ 
ditions. We can observe the phase transition from the 
low temperature regime, in which the correlation func¬ 
tion exhibits a power-law decay, to the high temperature 
regime, in which it displays an exponential decay. The 
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FIG. 1: (Color online) Spin correlation function C(r) for dif¬ 
ferent values of temperature. The upper (a) and lower (b) 
panels depict C(r) for the directions parallel and perpendic¬ 
ular to the driving velocity v, respectively. The values of 
temperature are T = 0.40, 0.60, 0.70, 0.80 and 0.90 from the 
top to the bottom of the panels. The other relevant parame¬ 
ters are ho = 1.0 and v = 0.5. The transition temperature is 
estimated as T c = 0.70±0.05. The symbols o, -|-, and x repre¬ 
sent C(r) for system sizes of 60 3 , 100 3 , and 150 3 , respectively. 
The error-bars are comparable with the size of the symbols. 
The panel (c) shows the temperature dependence of the ex¬ 
ponents oi|| and a± determined from the data corresponding 
to the system size of 60 3 . The error due to the fitting with 
a power function is displayed. The horizontal lines represent 
the values predicted from the spin-wave approximation. 


panel (c) represents the exponents as a function of tem¬ 
perature. The horizontal lines represent the theoretical 
prediction of the spin-wave approximation. The expo¬ 
nents ck|| and aj_ increase linearly at low temperatures, 
but they exhibit strong temperature dependence near the 
transition temperature. 

We next determine the transition temperature T c as a 
function of the strength of the random field ho and the 
driving velocity v by using the non-equilibrium relaxation 
method [2{|. With this method, we observe the relax¬ 
ation of the magnetization M(t) = V~ x f tf>(r)dr from 
the complete ordered state 0 1 (r*) = 1 and </> 2 (r) = 0. The 



FIG. 2: (Color online) Panel (a): Relaxation of the magne¬ 
tization M(t ) for different values of temperature. The other 
parameters are h o = 1.4 and v = 1.0. The values of temper¬ 
ature are T = 0.52, 0.55, 0.58, 0.65, and 0.75 from the up 
to the bottom. Panel (b): Scaling plot of the magnetization 
to Eq. (11511 with appropriately chosen t(T) and A. A = 0.08 
and r(0.75) = 1, r(0.65) = 4.0, r(0.58) = 18, r(0.55) = 
60, t(0.52) = 300. Panel (c): Relaxation time r as a function 
of temperature. The curve fitted to Eq. (1161) with T c = 0.43 
is shown. 


asymptotic behavior of the magnetization for T > T c is 
summarized as 


M(t) 


exp(-t/r(T)), (T > T c ), 
t~\ (T = T c ), 


(14) 


where t(T) is the relaxation time. In order to determine 
the relaxation time as a function of temperature for the 
disordered phase, we assume the following scaling form, 

M(t) = r(T)- A m(t/r(T)), (15) 

for T > T c . The relaxation of the magnetization M(t) 
and its scaling plot are displayed in the panels (a) and 
(b) of Fig. [2] The system size is 60 3 . 100 independent 
runs are performed for averaging. From the analogy to 
the KT transition of the two-dimensional pure XY model 
f30j , the correlation length is expected to diverge expo¬ 
nentially as £ ~ exp{a/s/T — T c ). Thus, we assume that 
the relaxation time diverges in the same way, 

r(T)=Bexp(—4=), ( 16 ) 

Fitting t(T) into Eq. (fl6l) with parameters A 1 B, and T c , 
we obtain the transition temperature. The best fitting is 
shown in the panel (c) of Fig. [2] 

Fig. [3] displays the schematic phase diagram with re¬ 
spect to the strength of the disorder ho, driving veloc¬ 
ity v. and temperature T. The QLRO phase appears in 
the large-D and low-T regime. In the region in which 
h 0 = 0, the LRO phase exists because the model is iden¬ 
tical to the three-dimensional pure XY model in the mov¬ 
ing frame. The small-u and high-T regime corresponds 
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FIG. 3: (Color online) Schematic phase diagram with respect 
to the disorder ho, driving velocity v, and temperature T. 
The parameter T c o denotes the transition temperature of the 
three-dimensional pure XY model. The abbreviation PM de¬ 
notes the paramagnetic or disordered phase. The inset depicts 
T c as a function of v. The solid lines serve as a visual guide. 
The squares (red) and circles (green) denote the values of T c 
for ho = 1.0 and ho = 1.2, respectively. 


to the disordered phase. The phase boundary at T = 0 is 
given by T/ig ~ Kv. It is noteworthy that the infinitesi¬ 
mally small random field breaks the LRO and leads to the 
QLRO. For an arbitrarily large value of ho, the QLRO 
is observed for sufficiently large values of v. We have 
checked that Eq. m is satisfied in the region of the 
QLRO phase displayed in the phase diagram Fig. [3] at 
low temperatures. This justifies the validity of the spin- 
wave approximation. 

In order to compare this transition with the KT transi¬ 
tion of the two-dimensional pure XY model, we calculate 
the specific heat as a function of temperature. We define 
the specific heat as c(T) = d(H) ss /dT , where (...} ss rep¬ 
resents the average with respect to the non-equilibrium 
steady state. Since the system is not in equilibrium, 
this specific heat is not related to the energy fluctua¬ 
tion. Fig. [4] shows the specific heat as a function of 
temperature for different values of the disorder strength 
ho- The upper arrows represent the transition tempera¬ 
ture T c as determined by the non-equilibrium relaxation 
method. The case that ho = 0 corresponds to the three- 
dimensional pure XY model. In the presence of a finite 
amount of disorder, the discontinuity of c(T) disappears 
and a smooth peak remains above the transition temper¬ 
ature. The position of the peak decreases with increase 
in the strength of the disorder. It is to be noted that the 
specific heat does not exhibit any singularity at the tran¬ 
sition temperature. The absence of a singularity and the 
existence of the smooth peak resemble the KT transition 
of the two-dimensional pure XY model Q ■ 



FIG. 4: (Color online) Specific heat c(T) as a function of 
temperature for system size of 60 3 . The driving velocity is 
v = 1.0. The crosses (red), squares (green), and circles (blue) 
denote c(T) values for ho = 0.0, 1.0, and 1.4, respectively. 
The upper arrows represent the transition temperature T c as 
determined by the non-equilibrium relaxation method. 



FIG. 5: (Color online) Transition temperature T c as a func¬ 
tion of v for the Heisenberg model (N = 3). The squares 
(red) and circles (green) denote T c values for ho = 0.8 and 
ho = 1.0, respectively. The solid lines serve as a visual guide. 


B. Heisenberg model 

Next, we consider the Heisenberg model (TV = 3). We 
found that the phase diagram of the Heisenberg model 
is qualitatively similar to that of the XY model. Fig. [5] 
shows the transition temperature determined by using 
the non-equilibrium relaxation method. We display the 
specific heat as a function of temperature in Fig. [G] In the 
presence of a finite amount of disorder, the discontinuity 
of c(T) disappears and a smooth peak remains above the 
transition temperature. The presence of the QLRO for 
N = 3 contrasts with the cases of the two-dimensional 
pure O (N) models, in which the KT transition does not 
exist for N = 3 pH] . 


V. CONCLUSIONS 

We investigated the non-equilibrium phase transition 
of three-dimensional O (N) models driven with a uniform 
velocity over a quenched random field. For the cases that 
N = 2 and 3 (XY and Heisenberg model, respectively), 
we show that QLRO appears in the strong driving regime. 
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FIG. 6: (Color online) Specific heat c(T) as a function of 
temperature for the Heisenberg model (N = 3). The driving 
velocity is v = 1.0. The crosses (red), squares (green), and 
circles (blue) denote c(T) values for ho = 0.0, 0.8, and 1.0, 
respectively. The upper arrows represent the transition tem¬ 
perature T c as determined by the non-equilibrium relaxation 
method. 


Furthermore, the transition from the QLRO phase to the 
disordered phase resembles the KT transition in the two- 
dimensional pure XY model. It is noteworthy that there 
exists a critical value N c such that QLRO does not ex¬ 
ist for N > N c . The renormalization group analysis is 
required to determine this value of N c . Moreover, it is 
also necessary to clarify the relation between this QRLO 
and the Bragg glass phase in the vortex lattice systems 
mm. We plan to investigate these theoretical aspects 
in our future studies. 

Finally, we remark on the topics related to the experi¬ 
mental realization of the QLRO under consideration. We 
consider nematic liquid crystals (NLCs) flowing in a ran¬ 
dom substrate such as an aerogel or a porous medium. 
Recently, the dynamics of liquid crystals confined to a 
complex geometry has attracted considerable attention 
due to not only fundamental research interest but also 
its industrial applications [HI, 33| . For NLCs in a ran¬ 
dom substrate, the random anchoring, which results from 
the complicated surface structure, significantly influences 
the ordering structure and thermodynamic properties 
[34l . S3. Here, we remark that certain factors concerning 
NLCs flowing in the random substrate are not included in 
our models. For example, we ignore the inhomogeneity of 
the velocity field. The inhomogeneity of the velocity field 
acts as an additional random perturbation for the direc¬ 
tors of the NLCs. Since the correlation of the spatial fluc¬ 
tuation of hydrodynamic velocity field exhibits power-law 
decay, this random perturbation is qualitatively different 
from the random anchoring whose correlation decays over 
a short distance. The consequences of this hydrodynamic 
effect and other factors excluded in our models should be 
investigated in future work. 
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Appendix A: Calculation of the spin-wave 
approximation 


We investigate the asymptotic behavior of the mean 
square relative displacement Eq. ©. First, we consider 
the case in which — r 2 is parallel to v. Eq. (JSJ) is then 
rewritten in terms of the polar coordinate 


= 2T 2 hl 


d q 


(27t) 3 j q j o 


<( 0 (n)- 0 (r 2 )) 2 ) 

27r p tc 

dc/) / df?g 2 sin 0 


x 


r 2 h 2 0 r°° 

27T 2 In 


d q 


1 — cos {g|rr — r 2 |cos$} 
F 2 K 2 q 4 + v 2 q 2 cos 2 9 
^ 1 - cos{g|ri - r 2 \u} 

_! T 2 K 2 q 2 + v 2 u 2 ’ 


(Al) 


where we have changed the variable u = cos 9. Since we 
are interested in the contribution from the small wave 
number regime q -C v/(TK), the range of the u-integral 
can be extended to (— 00 , 00 ). Thus, we have 


L 


v/(TK) x 

dg- 
q 


mn)-0{r 2 )) 2 ) 


1 — exp 



, r h 2 

2ttKv 

r 2 |)}. (A2) 


The integrand takes the maximal value at q = q m ~ 
y/v/(TK\ri — r 2 |), and it decays as q~ l for q q m . 
Therefore, we obtain 


mri)-0{r 2 )) 2 ) 

r/W dq 

2tt Kv J (hn q 


Th 2 




2irKv 


9m 


const 


const 


r^o , 

, „ login. ~r 2 \+ const. 
47 tKv 


(A3) 


This is one of Eqs. ©. 

We next consider the case in which ri — r 2 is perpen¬ 
dicular to v. Eq. © is then rewritten as 


mri)-9(r 2 )) 2 ) 

or 2 h 2 f d3q 1 ~ cos {g^l r i ~ r i\} 

^ (2tt ) 3 T 2 K 2 q 4 + v 2 q 2 x ’ 


(A4) 


where r\ — r 2 is assumed to be parallel to the 2 -axis. 
Since we are interested in the contribution from the small 
wave number regime q <C v/(TK) : the term containing 
q x in q 4 of the denominator can be ignored. Thus, we 
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have 


<(0(n)-0(r 2 )) 2 } 

f d 3 q 1 - cos{g z |ri - r 2 ]} 

q<v/(TK) ( 27r ) 3 r 2 K 2 (ql + g 2 ) 2 + u 2 g 2 ’ 


r/i-o f dq y d q z 1 - cos{g~|ri - r 2 |} 

ttA'z; J q<v/{ rK) 2t r 2 tt g 2 + g 2 

= r§_ r /(r * } dg l-Jp^ln-ral) 

7r/\'u 7 0 27 t g 

r/i 2 

~ ° log|ri - r 2 | + const, (A5) 

2nKv 


where Jq(x) is the Bessel function. This is the other one 
of Eqs. ©. 


Appendix B: Numerical calculation of the 
correlation function 


We explain the method of the numerical calculation 
of the correlation function. The correlation function is 
defined as 

C(t) = J DhP RF [h\ J D4>(j)(r) ■ 0(O)P st [</>; A], (Bl) 

where P st [0; h] denotes the steady state solution of the 
Fokker-Planck equation corresponding to Eq. © for a 
fixed random field h(r) and Prf[^] is the distribution 
function for h(r). 

We implement the numerical calculation in a moving 
frame with v. Thus, we define r' = r — vt and cj> (r r ) = 
4>(r). Eq. © is then rewritten as follows: 


dt 


6H[<j>’-,h t ] 




(B2) 


where h t {r’) = h(r' + vt) is a moving random field. 


The numerical calculation is implemented in a cubic 
simulation box, whose size is L. We numerically solve 
Eq. (1B2I) with the periodic and free boundary conditions 
for the directions perpendicular and parallel to the driv¬ 
ing velocity v, respectively. In what follows, we omit the 
prime symbols from Eq. (IB21) . The random field is contin¬ 
uously generated in the front boundary of the simulation 
box x = A, and it moves with velocity —v. h B (r±, t) de¬ 
notes the random field generated at the front boundary, 
where r± = (y,z) represents a coordinate of the front 
boundary. Then, the moving random field is written as 
h t { x, y, z) = h R {y , z,t—(L — x)/v). We denote the time 
dependent solution of the Fokker-Planck equation corre¬ 
sponding to Eq. (1B2I) for a realization of h R (r_i_,s) by 
P t CM [0; {/iB(s)}s<t]- Note that P t CM depends on all his¬ 
tory of the random field h R . From the equivalence of 
Eqs. © and ([Bl, Eq. dHU) is rewritten as follows: 

C( r ) = J Dh B P rf [M 

x J P00(r) • 0(O)P t CM [^;{h. B (s)} s < t ], (B3) 

for sufficiently large t, where Prf[Iib] is the distribu¬ 
tion function for the random field generated at the front 
boundary. 

If the simulation box is sufficiently large L 
where £ is the correlation length, we can assume the self¬ 
averaging property. Thus, the average with respect to 
the random held and the thermal noise can be replaced 
by the spatial average 

C(r) = y J dr 1 + r) ■ 4>(r’), (B4) 

where <j>{r) is a solution of the Langevin equation (IB2I) . 
We checked that the dependence of the correlation func¬ 
tions Eq. (1B41) on realizations of h B (r±,t) and r)(r,t) is 
small. We also take the time average of Eq. (1B41) at time 
intervals of L/v. 
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